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0^ ' We consider the 2D Hubbard model on the honeycomb lattice, as a model for a single layer 

, graphene sheet in the presence of screened Coulomb interactions. At half filling and weak enough 

' coupling, we compute the free energy, the ground state energy and we construct the correlation 

functions up to zero temperature in terms of convergent series; analiticity is proved by making use 
f— I , of constructive fermionic renormalization group methods. We show that the interaction produces a 

■ modification of the Fermi velocity and of the wave function renormalization without changing the 
asymptotic infrared properties of the model with respect to the unperturbed non-interacting case; 

, this rules out the possibility of superconducting or magnetic instabilities in the ground state. We 

■ also prove that the correlations verify a Ward Identity similar to the one for massless Dirac fermions, 
up to asymptotically negligible corrections and a renormalization of the charge velocity. 

S: 

I. INTRODUCTION 

, The recent experimental realization of a monocrystalline graphitic film, known as graphene 

[2d|, revived the interest in the low temperature physics of two-dimensional electron systems 
on the honeycomb lattice, which is the typical underlying structure displayed by single-layer 
graphene sheets. Graphene is quite different from most conventional quasi-two dimensional 
' electron gases, because of the peculiar quasi-particles dispersion relation, which closely resem- 

^ , bles the one of massless Dirac fermions in 2 + 1 dimensions. This was already pointed out in 

[25I and further exploited in , where the analogy between graphene and 2 + 1-dimensional 
quantum electrodynamics (QED) was made explicit, and used to predict a condensed-matter 
analogue of the axial anomaly in QED. From this point of view, graphene can be considered as 
a sort of testing bench to investigate the properties of infrared QED in 2 -I- 1 dimensions. Re- 
cently, the experimental observation of graphene greatly enhanced the study of the anomalous 
' effects induced by the pseudo-relativistic dispersion relation of its quasi particles, see Q for an 

, up-to-date description of the state of art. Among the most unusual and exciting phenomena 

displayed by graphene, and already experimentally observed, let us mention the anomalous 
integer quantum Hall effect and the insensitivity to localization effects generated by disorder. 
It is reasonable to guess that the unique properties of graphene will have in the next few years 
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\ several important applications in condensed matter and in nano-technologies. 

The main reason behind these anomalous effects lies in the geometry of the Fermi surface, 
which at half filling is not given by a curve, as in usual 2D Fermi systems, but is completely 
degenerate: it consists of two isolated points^ as in one dimensional Fermi systems. From a 
theoretical point of view, this fact completely changes the infrared scaling properties of the 
propagator. It has been pointed out, see for instance |llj] and references therein, that, in the 
case of short-range electron-electron interactions, all the operators with four or more fermionic 
fields are irrelevant in a Renormalization Group (RG) sense; this suggests that the interaction 
should not affect too much the asymptotic behavior of the model, at least at small coupling. 
It should be remarked however that such RG analyses were performed only at a perturbative 
level, without any control on the convergence of the expansion, and directly in the relativistic 
approximation, consisting in replacing the actual dispersion relation by its linear approximation 
around the singularity; such approximation implies in particular the validity of a continuous 
Lorentz U{1) symmetry that is not present in the original model. 

Aim of this paper is to present the first rigorous construction of the low temperature and 
ground state properties of the 2D Hubbard model on the honeycomb lattice with weak local 
interactions; this is achieved by rewriting the correlation functions in terms of resummed series. 
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convergent uniformly in the temperature up to zero temperature, as we prove by making use 
of constructive fermionic renormalization group. We show that indeed the interaction does 
not change the asymptotic infrared properties of the model with respect to the unperturbed 
non-interacting case, but it produces a renormalization of the Fermi velocity and of the wave 
function (note that no renormalization of the Fermi surface would be present in the relativistic 
approximation). Our result rules out the presence of superconducting or magnetic instabilities 
at weak coupling; this is in striking contrast with the Hubbard model on the square lattice, 
where quantum instabilities (corresponding to the magnetic or superconducting long range order 
that are presumably present in the ground state) prevent the convergence of the perturbative 
expansion in U for low enough temperatures. We also prove that indeed the 2D Hubbard 
model on a honeycomb lattice is asymptotically described by a QED2+1 in the presence of an 
ultraviolet cutoff, massive "photons" and massless electrons; however the bare parameters of 
the QED theory must be carefully chosen to include lattice effects. 



II. THE MODEL AND THE MAIN RESULTS 



A. The model 



The grandcanonical Hamiltonian of the 2D Hubbard model on the honeycomb lattice at half 
filling in second quantized form is given by: 

iJA = - V V fa+ fo" - +6+ - aZ ) + (2.1) 



1 = 1,2,3 



where: 



1. A is a periodic triangular lattice, defined as A = M/LM, where L G N and B is the 
triangular lattice with basis ai = ^(3, V3), 0,2 — 5 (3, — %/3). 

2. The vectors bi are defined as 

Ji = (l,0), J2 = i(-l,V3), l3 = i(-l,-V3) . (2.2) 



3. ai^ are creation or annihilation fermionic operators with spin index a —\\, and site index 
af £ A, satisfying periodic boundary conditions in x 

4. 6^ -, are creation or annihilation fermionic operators with spin index cr =ti and site 
index x + 5; G A + (5i, satisfying periodic boundary conditions in x. 

5. f7 is the strength of the on-site density-density interaction; it can be either positive or 
negative. 

Note that the Hamiltonian (j2.ip is hole-particle symmetric, i.e., it is invariant under 
the exchange ai < — *-a5 , 6* - < — > — ^ . This invariance implies in particu- 

lar that, if we define the average density of the system to be p — (2|A|)^^ (A^)^ ^, 
with — „(a% „aZ -f 6^ -. 6Z ? ) the total particle number operator and 

(•)^ = Tr{e~'^^'^-}/Tr{e~'''^'^} the average with respect to the (grandcanonical) Gibbs 
measure at inverse temperature /3, one has p = 1, for any |A| and any fi. We also recall that 
a theorem [l^ guarantees that at half filling the ground state of (|2.ip is unique and its total 
spin is equal to zero. 



Our goal is to characterize the low and zero temperature properties of the system described 
by (j2.ip . by computing thermodynamic functions (e.g., specific free energy and specific ground 
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state energy) and a complete set of correlations at low or zero temperatures. To this purpose 
it is convenient to introduce the notions of specific free energy 

/^(t/) = -i lim |ArilogTr{e-^^n, (2.3) 

P |A|^oo 

of specific ground state energy e[U) — liuip^oo fpiU), and of Schwingcr functions, defined as 
follows. 

Let us introduce the two component fermionic operators „ — {(^3 ) ^^'^ 

us write Vpi , = ai and ^'i „ — bt , ■ We shall also consider the operators — 
gHxoxjf^ ^g-Hxo X = {xo,x) and xq e [0,/?], for some /3 > 0; we shall call xq the time 

variable. We shall write ^E"^^ ^ = and ^'^ j 2 = ^x+Si o-' ^^^^ ^1 = (0: ^i)- We define 

5f'^(xi,ei,ai,pi;...;x„,e„,a„,p„) - {T{^'^\^,^^p^ ■ ■ - ^'^l^^^^pj) (2.4) 

where: x; G [0, /?] x A, (Ti =ti, £i = ±, Pi = 1,2 and T is the operator of fermionic time 
ordering, acting on a product of fermionic fields as: 

nn\,a,.p, ■ ■ ■ *xi.„,pj = (-i)'^c:'<?>,..„)-p.(i) ■ ■ ■ *x:t'„^.(„),p.(„) (2.5) 

where tt is a permutation of {1, . . . ,n}, chosen in such a way that 0:7^(1)0 ^ ■ ■ ■ ^ 2;^(„)o, s-i^d 
(— 1)'^ is its sign. [If some of the time coordinates are equal each other, the arbitrariness of 
the definition is solved by ordering each set of operators with the same time coordinate so that 
creation operators precede the annihilation operators.] 

Taking the limit A — > 00 in (j2.4p we get the finite temperature n-point Schwinger functions, 
denoted by ^^(xi, ei, cti, pi; . . . ; x„, e„, (T„, p„), which describe the properties of the infinite 
volume system at finite temperature. Taking the /3 — > 00 limit of the finite temperature 
Schwinger functions, we get the zero temperature Schwinger functions, simply denoted by 
5'„(xi,ei, (Ti, pi; . . . ; x„, £„, (7„, p„), which describe the properties of the ground state of (j2.ip in 
the thermodynamic limit (note that in this case, by the uniqueness of the ground state proved 
in [l6l |. the infinite volume and zero temperature limits commute). 



B. The non interacting case 



In the non-interacting case U = the Schwinger functions of any order n can be exactly 
computed as linear combinations of products of two-point Schwinger functions, via the well- 
known Wick rule. The two-point Schwinger function itself, also called the free propagator, for 
X 7^ y and x — y 7^ (±/3, 0), is equal to (see Appendix [Al for details): 



So'^{x-y)p^p' = 5^'^(x,CT, -,p;y,cr, +,p') 



U=Q 



M — irv-> 



M~^oo(3\A\^^^^^kl + \v{k)\^\~v{k) tko 
where: 

1. 7\f e N, k = (fco, k) and V^^l ^Vpx Vl; 

2. - {ko = f (no + i) : no = -Af, . . . , M ~ 1} and Vl = {k = + : < 
n-i,n2 < L~l}, where bi = ^(1, V3), ^2 ~ if (1; ^V3) are a basis of the dual lattice A*; 

3. v{k) = ^3^^ e*'=(^'~^i) = 1 + 2e-*3/2fei cos^fca; its modulus |'i;(fc)| is the dispersion 
relation, given by 



= ^J (1 + 2cos(3fci/2) cos(\/3fc2/2))^ + 4sin2(3fci/2) cos2(V3fc2/2) . (2.7) 
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Atx = yorx — y = (±/3,0), the free propagator has a jump discontinuity, see discussion 
at the end of Appendix [X\ Note that S'q'^(x) is antiperiodic in xq, i.e. Sq''^{xo + P,x) = 
—Sq'^{xo,x), and that its Fourier transform 6*0 '^(k) is well-defined for any k e 2?/3,l, even in 
the thermodynamic limit L 00, since j/col > ^- We shall refer to this last property by saying 
that the inverse temperature /3 acts as an infrared cutoff for our theory. 

If we take (3, L 00, the limiting propagator 5*0 (k) becomes singular at {ko — 0}x{fc = jjj}, 
where 

— + /27r 27r . , , 

are the Fermi points (also called Dirac points, for an analogy with massive QED2+1 that will 
become clearer below). Note that the asymptotic behavior of v{k) close to the Fermi points is 
given by v{jy^ + k') ~ ^{ik'i ± k'2). In particular, if w = ±, the Fourier transform of the 2-point 
Schwinger function close to the Fermi point jJf, can be rewritten in the form: 

Qji,, I -*^o -vf{~ik[+ujk'^)+r^{k')\ 

where Zq = 1 is the free wave Junction renormalization and v^p' = 3/2 is the free Fermi velocity. 
Moreover, < Clfc'p, for small values of k' and for some positive constant C. 



C. The interacting case 

We are now interested in what happens by adding a local interaction. In the case 
[/ 7^ 0, the Schwinger functions are not exactly computable anymore. It is well-known 
that they can be written as formal power series in [/, constructed in terms of Feynmann 
diagrams, using as free propagator the function 5o(x) in (j2.6p . Our main result consists in a 
proof of convergence of this formal expansion for U small enough, after the implementation of 
suitable resummations of the original power series. Our main result can be described as follows. 

Theorem 1. Let us consider the 2D Hubbard model on the honeycomb lattice at half filling, 

defined by i2. There exist a constant Uq > such that, if \U\ < Uq, the specific free energy 
fpiU) and the finite temperature Schwinger functions are analytic functions of U , uniformly 
in [3 as (3 —f 00, and so are the specific ground state energy e{U) and the zero temperature 
Schwinger functions. The Fourier transform of the zero temperature two point Schwinger func- 

def 

tion S{x)p^pi = 5'2(x, (7, — , p; 0, (T, +, p'), denoted by 5'(k), is singular only at the Fermi points 
k = = (OyPp), see i2.8\) . and, close to the singularities, ifoj — zL, it can be written as 

+ "'5) + , (2.10) 

with k' = {kQ,k'), and with Z and vp two real constants such that 

Z = 1 + a[/2 + 0(C/3) , VF = ^+bU^ + 0{U^) 

where a and b are non-vanishing constants. Moreover the matrix -R(k') 
||^(k')|| ^ Cjk'l'' for some constants C,^ > and for |k'| small enough. 

Remarks. 

1) Theorem 1 says that the location of the singularity does not change in the presence of 
interaction; on the contrary, the wave function renormalization and Fermi velocity are modified 
by the interaction. Note also that, in the presence of the interaction, the Fermi velocity remains 
the same in the two coordinate direction even though the model does not display 90° discrete 
rotational symmetry, but rather a 120° rotational symmetry. 
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2) The resulting theory is not quasi-free: the Wick rule is not valid anymore in the presence of 
interactions. However, the long distance asymptotics of the higher order Schwinger functions 
can be estimated by the same methods used to prove Theorem 1, and it is the same suggested 
by the Wick rule. 

3) The fact that the interacting correlations decay as in the non-interacting case implies in 
particular the absence of long range order at zero temperature, e.g., the absence of Neel order 
in the ground state at weak coupling. In fact, as a corollary of our construction, we find: 

lim {Ss ■ Sg) A < C—^ , (2.12) 

where, if x £ A, the spin operator 5*^ is defined as: Sg = ^s.'^'^s with ai, i = 1,2,3, the 
Pauli matrices; similarly, if x G A + 5i, = J^a- ■'^^s ■■ ^ot^ that it is well known that the 
ground state has zero total spin (iGj , however existence of Neel order was neither proven nor 
ruled out by the results in [ig . 

4) Similarly to what remarked in the previous item, one can exclude the existence of super- 
conducting long range order: the Cooper pairs correlations decay to zero at infinity at least as 
fast as the spin-spin correlations in (|2.12p . 

5) Our analysis can be extended in a straightforward way to the case of exponentially decaying 
interactions (instead of local interactions). However, if the decay is slower, the result may 
change. In particular, in the presence of 3D Coulomb interactions, the electron-electron 
interaction becomes marginal (instead of irrelevant), in a renormalization group sense p^ . 

6) Previous analyses of the Hubbard model on the honeycomb lattice were performed only at 
a perturbative level, without any control on the convergence of the weak coupling expansion, 
and directly in the Quantum Field Theory approximation, consisting in the replacement of 
5'o(k) by its linear approximation around the Fermi points, sec for instance [ll| and references 
therein. 

7) In Appendix [C] we prove that indeed the 2D Hubbard model on a honeycomb lattice is 
asymptotically described by a QED2+1 model in the presence of an ultraviolet cutoff, massive 
"photons" and massless electrons, provided that the bare parameters of the QED2+1 are 
carefully chosen to include lattice effects. As a result, the correlations asymptotically verify a 
modified Ward Identity (WI) related to an approximate local U{1) Lorentz symmetry: note, 
however, that the renormalized charge velocity appearing in the modified WI for graphene 
explicitly breaks rotational invariance, contrary to what happens to the Fermi velocity, or to 
the charge velocity of a pure relativistic QED model. 

The proof of the Theorem is based on constructive fermionic Renormalization Group (RG) 
methods, see [2l.[l8l[22j for extensive reviews. It is worth remarking that the result summarized 
in Theorem 1 is one of the few rigorous construction of the ground state properties (including 
correlations) of a weak coupling 2D Hubbard model. The only other example we are aware of is 
the Fermi liquid construction in [3] , applicable to cases of weakly interacting 2D Fermi systems 
with a highly asymmetric interacting Fermi surface. Related results include the construction of 
the state at temperatures larger than a BCS-like critical temperature 0,01) or the computation 
of the first contribution to the ground state energy in a weak coupling limit (lol . [TtI . [23| . 

The rest of the paper will be devoted to the proof of Theorem 1. In Sec. IIII AI we review 
the Grassmann integral representation for the free energy and the Schwinger functions. In 
Sec lIII Bj we start to describe the integration procedure leading to the computation of the free 
energy, and in particular we describe how to integrate out the ultraviolet degrees of freedom. In 
Sec lIII C] we complete the proof of convergence of the series for the free energy and the ground 
state energy. In Sec lIII Dl we describe the proof of convergence for the series for the Schwinger 
functions, with particular emphasis to the case of the two-point Schwinger function. In the 
Appendices we provide further details concerning the non-interacting theory, the ultraviolet 
integration and the equivalence (as far as the long distance behavior is concerned) between the 
Hubbard model and a massive QED theory in 2+1 dimensions. 
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III. RENORMALIZATION GROUP ANALYSIS 



A. Grassmann Integration 

It is well-known that the usual formal power series in U for the partition function and for 
the Schwinger functions of model (j2.ip can be equivalently rewritten in terms of Grassmann 
functional integrals, defined as follows. 

We consider the Grassmann algebra generated by the Grassmannian variables 
{^t<yjkevir^''^ ^ Grassmann integration / [Ukev,,r.Ua^n di^+^^d-^^^^^J defined 
as the linear operator on the Grassmann algebra such that, given a monomial (3(^^,4'^) 
in the variables *I'kcrp, its action on is except in the case (5(4'^, ^P^) = 

Yluev^ lUcZ^U a p^i. a ^ pcrmutatiou of the variables. In this case the value 

of the integral is determined, by using the anticommuting properties of the variables, by the 
condition 

/[ n Ud^LjK..] n UK^A.=^ (3-13) 

Defining the free propagator matrix as 

au-i^'k (3.14) 



-v{k) —iko 
and the "Gaussian integration" P{d'i') as 



• exp{ - iP\A\r' K^Jk'K^^.} , (3.15) 



'^=11 



it turns out that 



^(rf*)*k„.,,p,*^.,..,p. =/3|A|5....<5k,,k.[5kJp^,p^ , (3.16) 
so that, if X - y ^ /3Z X {0}, 

J™ -Jji E e-'^'^-y^ffk = lim f Pm^Z.A = Soi^-y) , (3.17) 

M^oc fj\A\ M^ooJ 

kfciJ;3,t 

where S'o(x — y) was defined in p.5p and the Grassmann fields ^'x a ^-''^ defined by 

"^i^^p^WXl E ^^^''^t.P' xeA^,MxA, (3.18) 

with A/3,A/ = {m/3/M : m ~ —M,...,M — 1}. Let us now consider the function on the 
Grassmann algebra 



v{^) ^ uJ2 f d^^ltpK^.p^^tupK^p 

^ — 1 O 



p=l,2 

u 



f«|A|)3 E E *^p,T,P*k.T.p*^'+p,i.p*k',i,p , (3.19) 

p=l,2 k,k',p 
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where, in the first hne, the symbol J dx must be interpreted as 

^-=4 E' (3-20) 

and, in the second hne, the sums over k, k' run over the set I?/3.l, while the sums over p run over 
the set 27r/3~^Z x Vl (p is the transferred momentum). Note that the integral J P((i4')e^^^*^ 
is well defined for any U; it is indeed a polynomial in U, of degree depending on M and L. 
Standard arguments show that, if there exists the limit of / P((i5')e^^'*-' as AI oo, then the 
normalized partition function can be written as 

e^/3|A|F,,. dg- Tr[e-^-^^] ^ /p(dvl/)e-^(*) (3.21) 

where Hq is equal to p.l|) with U — 0. A possible way to prove (13.211) is to compare the 
perturbation theory obtained by expanding in powers of U via Trotter's product formula the 
trace Tr{e-^^'^}/Tr{e-'^^«} with the one obtained by expanding in U the Grassmann func- 
tional integral, and then show that they are the same, order by order, see [ij. This proof also 
shows that the correct choice of the interaction (|3.19p expressed in Grassmann variables does 
not include terms bilinear in the fields, contrary to the interaction in second quantized form, 
see (|2.ip : in fact, with this choice, in both perturbative expansions the "tadpoles" are exactly 
vanishing, as required by the condition that the system is at half filling. 

Similarly, the Schwinger functions at distinct space-time points, defined in (|2.4p . can be 
computed as 

/P(d«')e-^(*)^'Ji ^ „ ••■*x" <r „ 
S{xi,(Ti,ei,pi; . . . ;x„,cr„, £„,/)„) = hm^ / P(d^)e-^('^) ^^^^^ ■ (3.22) 

Note that the limit x.^ — x^ and the limit M — > oo do not commute in general. 

In the following we shall study the functional integrals by introducing suitable expansions 
where the value of M plays no essential role and we shall indeed be able to control such 
expansions uniformly in M , if U is small enough and that the limit M oo can be taken 
in the resulting expressions for the free energy and the Schwinger functions. For this reason, 
from now on we shall not stress anymore the dependence on M, unless for the cases where the 
presence of a finite M is relevant, e.g., for the analysis of the ultraviolet integration described 
in Appendix [XI 

It is important to note that both the Gaussian integration P{d'^) and the interaction V{'i>) 
are invariant under the action of a number of remarkable symmetry transformations, which 
will be preserved by the subsequent iterative integration procedure and will guarantee the 
vanishing of some running coupling constants (see below for details). Let us collect in the 
following lemma all the symmetry properties we will need in the following. 

Lemma 1. For any choice of M, f3, A, both the quadratic Grassmann measure P{d'^) defined 
in I13.15\) and the quartic Grassmann interaction V{^) defined in h3.19[) are invariant under 
the following transformations: 

(1) spin exchange : *k,o-,p^^'k,-^,p/ 

(2) global U{1): '^l, „ p —> e''^°''''^l,,„,pj ao- G M independent o/ k; 

(3) spinSO{2): ( *^^T.p ] ^ ^ [ *k,T,P ] yj^th r ^ ( cos0 sin0\ q ^ j 

— — - \n,i,J \n,i,J y-sme cose] 

independent of k; 

(4) discrete spatial rotations: - ^ e^F^^^-'^-^Olp-i)^* , with Tix"^^^ R2^/3x; 
note that in real space this transformation simply reads ^ ^^-^ ^ and 

(5) complex conjugation: (j p ^ up' c ^ c* , where c is a generic constant appearing in 
P{d^) and/or in F(*); 

(6.a) horizontal reflections : ^fk„k„k,).a.i^'^tko.~k,.k,).^.2' 
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(6.b) vertical reflections : ^koMM).^,P ^ ^tk»M.-k^),<y,P' 

(7) particle-hole: i#7 

(/co,fi:),o',p [kQ, — k).a,p 

(8) inversion : - —f i(-l)Pii^ 

Proof. A moment's thought shows that the invariance of V(\E') under the above symmetries 
is obvious, and so is the invariance of P{d^) under (l)-(2)-(3). Let us then prove the invariance 
oi P{d^) under (4)-(5)-(6.a)-(6.b)-(7)-(8). More precisely, let us consider the term 

E*t..-5k'*k..,. - (3-23) 

k 

E ^k,..i^o^k..4 - E ik)i>^,.a - E ^k,.,2^(fc)^k,.4 - * E *i.,2fco^k,.,2 



in (|3.15p . and let us prove its invariance under the transformations (4)-(5)-(6.a)-(6.b)-(7)-(8). 

Under the transformation (4), the first and fourth term in the second line of p.23p are 
obviously invariant, while the sum of the second and third is changed into 

^ L (ko.Tik).a,l ^ ' (ka.Txk).a.1 {ka .Txh) ,a .1 ^ ' (ko,Tik),cr,l 

k 

= ^ E [*J,.4^*(rr'^)e^''^'^"'^^*k,.,2 + *^,.,2e""'^'^"'=^«(rr'^)^k,.,i] ■ (3-24) 

k 

Using that v{T^^k) = e'''^^^-^'''>v(k), as it follows by the definition v{k) = Ei=i,2.3 

we find that the last line of (|3.24p is equal to the sum of the second and third term in p.23p . 

as desired. 

The invariance of (|3.23p under the transformation (5) is very simple, if one notes that v(^k) = 
v*{k), as it follows by the definition of v{k). 

Under the transformation (6. a), the sum of the first and fourth term in the second line of 
p.23p is obviously invariant, while the sum of the second and third is changed into 

" E ^(feo:-*:i,fc2),o-,2^ (^)^(feo,-fei,fc2),o-4 ^ E *^(feo:-*:i,fc2),o-4^'- -'^(feo, -'£1,^2), £7,2 ~ 
k k 

= -E*t^X((-^i'^2))«'k,.,i-E*t^.i<(-^i'fc2))«'k,.,2 ■ (3-25) 

k k 

Noting that v{{—ki, k2)) — f*(k), one sees that this is the same as the sum of the second and 
third term in (I3.23p . as desired. 

Similarly, noting that —k2)) — w(k), one finds that p.23p is invariant under the trans- 

formation (6.b). 

Under the transformation (7), the sum of the first and fourth term in (|3.23p is obviously 
invariant, while the sum of the second and third term is changed into 

+ y'I'" - v*{k)i!+ - +y*" v(k)^+ - 

(ka-k),a.l ^ ' (ko-k),a.2 Z-^ {ko.-k).cr,2 ^ ' (ko.-k).cr,l 
k k 

- - E (-fc)*k...i - E *J..,i^(-^)*k..,2 ■ (3-26) 

k k 

Using, again, that v{-~k) — v*{k), we see that the latter sum is the same as the sum of the 
second and third term in (|3.23p . as desired. 

Finally, under the transformation (8), all the terms in the right hand side of (|3.23p are 
separately invariant, and the proof of Lemma 1 is concluded. ■ 

B. Free energy: The ultraviolet integration 

We start by studying the partition function 

Ep^L = e-'^l^l^'^'^ = / P(d*)e^^(*) . (3.27) 
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Note that our lattice model has an intrinsic ultraviolet cut-off in the k variables, while the 
/cq variable is unbounded. A preliminary step to our infrared analysis is the integration of 
the ultraviolet degrees of freedom corresponding to the large values of ko- We proceed in the 
following way. We decompose the free propagator into a sum of two propagators supported 
in the regions of fco "large" and "small" , respectively. The regions of fco large and small are 
defined in terms of a smooth support function xoit) which is 1 for t < ao and for t > ao7, 

7 > 1; ao is chosen so that the support of xo (^\/ kg + \k — p^\^ ^ and xo (\/ + 1^ ^ PfP ) 

disjoint (here | • | is the euclidean norm over M^/A*). In order for this condition to be satisfied 
it is enough that 2ao7 < 47r/(3\/3); in the following, for reasons that will become clearer later 
we shall assume the slightly more restrictive condition 2ao7 < 47r/3 — 47r/(3V3). We define 



are 



/„.„.(k) = 1 - Xo (V^o + 1^ - P^\^ ) - Xo(V^o + \k-Pp\^) (3.28) 
and /i.r. (k) = 1 — /.tj.i,.(k), so that we can rewrite as: 

5k = /n... (k)gk + (k)5k 5<"-"-Hk) + (k) . (3.29) 

We now introduce two independent set of Grassmann fields {^^i^'^p^} and {^'k^p^}, with 
keVpL,o- =ti, p = 1,2, and the Gaussian integrations P(d*(" " )) and P(d*(* '')) defined 

by 

J p(rf*(--))*Lr.;\;,*Lr.i:, =/3iAi5.,,.,5k. , 

I Pid^^'-'^-W^ZL^t^JL = /3|A|<5.„.,<5k„k.5'"-)(ki)p„p. . (3.30) 

Similarly to P(d5'), the Gaussian integrations F(d5'(" "-)), F(d4'(' '"-^) also admit an expHcit 
representation analog ous to ([SH]) . with 5k replaced by g(" " )(k) or ^('■'^ ^(k) and the sum over 
k restricted to the values in the support of /u.t>.(k) or (k), respectively. It easy to verify 
that the ultraviolet propagator g("-''-)(x - y) = (/3|A|)-i J2i,eVff ^ e-*''(^-y)g("-"-)(k) satisfies 

|g'"-^-^(^-y)l< i + l^'!y|A. - (3-31) 

The definition of Grassmann integration implies the following identity ("addition principle"): 
P(d*)e-^W ^ J P(d*(* '' )) J p(d*(" -))e-^(**"' (3.32) 
so that we can rewrite the partition function as 

S,,, = e-'^l^l^-^^ /p(dvI;(-))exp{^l£:,^„ (-F(*(-) + .).„)} ^ 

n>l 

= e-^'l^l^" J P(dvI/(^-'^-))e-^(*'" ') , (3.33) 

where the truncated expectation fj^, is defined, given any polynomial Vi(^^"'^'') with coeffi- 
cients depending on as 

£L.iVi{-y,n) - ^log / P(dM'(-))e^^^(*'""') , ^ (3.34) 



A=0 



and V is fixed by the condition V(0) = 0. It can be shown (see discussion after (|3.37|) and 
Appendix [B)) that V can be written as 



C30 n 

v(*) = E(/?|Air" E E E [n 

n=l o-i,...,(T„ = Ti Pi.---'P2n = l,2 ki,...,k2„ J = l 



K2j-1 ,p2j-l K2j ,(Tj ,p2j 



•W^2n,p(kl, . . . ,k2„-l) (5(^(k2,_l - k2j)) , (3.35) 
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where p = (pi, . . . , p2ri) and we used the notation 

<5(k) = 5(fc)<5(fco) , '5(fc)-|A| ' ^(fco) - /34o,o , (3.36) 

ni,n2SZ 

with 61 , 62 a basis of A* . The possibihty of representing V in the form (|3.35p , with the ker- 
nels W2n,p independent of the spin indices Ui, follows from the symmetries listed in Lemma 1 

and from the remark that P((i4''-" " -') and P((i5'(* ''-^) are separately invariant under the same 
symmetries. 

The constant Fq in (|3.33|) and the kernels W2n,p in (|3.35p are given by power series in [/, 
convergent under the condition \U\ < Uq, for Uq small enough; after Fourier transform, the 
x-space counterparts of the kernels W2n,p satisfy the following bounds: 

/'dxi...dx2„[ n |x,-x,r-^]|l¥2„,p(xi,...,x2„)| </3|A|C;;|;7r^'^{i'"-i> , (3.37) 



<i<j<2n 



l<i<j<2n 

for some constant Cm > 0, where m = '^-^ 

The proof of convergence of the power series defining Fo and M^2n,p, as well as the proof of 
the bounds (|3.37p . uses the decay property p.3ip combined with standard fermionic cluster 
expansion methods. The proof is much simpler than the infrared integration that we shall 
study below, and is based on similar ideas; see Appendix [B] for a proof. Note that the decay 
(|3.3ip suggests the possibility of a single scale integration of the ultraviolet degrees of freedom. 
However, the discontinuity of the propagator sX xq — yo — implies that does not admit 
a Gram representation j2l| and this fact prevents the direct implementation of a single step 
fermionic cluster expansion (see next section for the notion of Gram determinant and for a 
description of the use of Gram determinants in the infrared multiscale integration) . A possible 
way out of this problem is to decompose the ultraviolet propagator as a sum of propagators, 
each admitting a Gram representation, and to perform a simple multiscale analysis of the 
ultraviolet problem. This strategy was described many times before in the literature, see 
for instance [1, 0, H, H, [3; for completeness, it will be presented in a self-contained form 
in Appendix |BJ Recently, a different proof based on a single scale integration step and us- 
ing improved bounds on determinants associated to "chronological products" was proposed [2T| . 

It is important for the incoming discussion to note that the symmetries listed in Lemma 
1 also imply some non trivial invariance properties of the kernels. We will be particularly 
interested in the invariance properties of the quadratic part V[^2,(pi,p2)('^)' which will be used 
below to show that the structure of the quadratic part of the new effective interaction has 
the same symmetries as the free integration. The crucial properties that we will need are the 
following. 

^ Lemma 2. Let VF,,(k) = VFs.ci.dW, Wbb{k) = VK2,(2,2) (k), WabCk) = W^2,(i,2)(k) and 
Wf,o(k) — W2 (2 i)(k). Then the following properties are valid: 

(l) Waai'k) = V;,b(k) and Wabi'k) = VFftKk),- 

(ii) as 13^ 00, for u = ±, W^aa(0,P^) = VF,b(0,p^) = 0; 
(Hi) as /3, |A| — > 00, for uj — ±, 

d^Waa{Q,PF) = , Re{dk,Waa{^),PF)] = , 9fco VFab (0, ) = , (3.38) 

Re{dk,Wab{QJF)] = lm{dk,Wab{Q.PF)] = , idk,Wab{Q,PF) = U0dk,Wab{Q,PF) ■ 



Remarks. 

1) For simplicity, the properties (ii) and (iii) are spelled out only in the zero temperature limit 
and in the thermodynamic limit; however, as it will be clear from the proof, those properties 
all have a finite temperature/ volume counterpart. 

2) Lemma 2 implies that in the vicinity of the Fermi points the kernel W2^(p^p')(k) can be 
rewritten in the form 
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for some real constants zo,5o, modulo higher order terms in {ko,k'). Therefore, it is apparent 
that its structure is the same as the one of S'o(k), modulo higher order terms in (fco, k'). 

Proof. As remarked after (j3.36|) . P{d^^'^ '"-^) and are separately invariant under 

the symmetry properties listed in Lemma 1. Therefore V('I') is also invariant under the same 
symmetries, and so is the quadratic part of V(\l/), that is 

k,p 

+ ^£-;i^*S;,2w^^^(k)] . (3.40) 

Recall that, as assumed in the lines preceding (|3.28[) . the support of 'I''^*''-) consists of two disjoint 
regions around andp^, respectively; in particular, we assumed that 2aoj < in / 3— An / (3^/3) . 
Under this condition, it is easy to realize that if both k and p+k belong to the support of 
then IpI < An/3. As a consequence, in (13.371) . the only non zero contributions correspond to 
the terms with p = (in fact, if p is 7^ and belongs to the support of 6{p), then |p| > 47r/3, 
which means that either k or k + p is outside the support of '!'(*•''■), and the corresponding term 
in the sum is identically zero). This means that the sum 



<T,k 



+*L:;;^'i'E;i"w^.a(k) + ^^:^^n:^^wM\ . (3.41) 

is invariant under the symmetries (l)-(7) listed in Lemma 1. 
Invariance under symmetry (4) implies that: 

Waa{ka,k) = Waa{ko,T-'k) , Wtb{ko, k) = Wtb{ko, k) , (3.42) 
Wab{ko,k) = e'^(''~'''>Wab{ko,T^'k) , Wba{ko,k) = e-'^(''-'''>Wab{ko,Tr'k) ; 

invariance under (5) implies that: 

Waa(k) = Waai-'k.y , Wbhik) = T^,,6(-k)* , (3.43) 
VKafc(k) = Wabi-'k)* , WbaCk) = W^6a("k)* ; 

invariance under (6. a) implies that: 

Waa{ko,ki,k2) = Wbb{k(),-ki,k2) , Wab{ko,ki,k2) = Wbaiko,-ki,k2) ; (3.44) 

invariance under (6.b) implies that: 

Waa{ko,kl,k2) = Waa{ko,ki,-k2) , Wfef, (fco , /ci , fc2 ) = W^66 (fco , fci , -fc2) , (3.45) 
Wabiko, fcl, fc2) = Wab(fco, fcl, -^2) , Wbaiko, fcl, fc2) = Wha{ko, fcl, -fc2) ; 

invariance under (7) implies that: 

Waaiko, k) = Waa{ko, -k) , W,,6(fco, fc) = Wbb{ko. "fc) , (3.46) 
Wab{kQ,k) = Wba{kQ,-k) ; 

Finally, invariance under (8) implies that: 

Waaiko, k) = -VK„a(-fco, fc) , Wbb{ko,k) = ~Wbb{-ko, fc) , (3.47) 
Wab{ko,k) = Wabi-ko, fc) , Wba{ko, k) = Wba{-ko, fc) ; 

Now, combining the first of (I3.44|) . the second of ()3.45p and the second of p.46p . we find that 
WaQ(k) = T/Ff,6(k). Combining the third of dSlSl), the third of (P^ISD and the last of (EUD, we 
find that Wabi'k) = W4a(k)*. This concludes the proof of item (i). 
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The first of (|3.47p implies that, as /3 ^ oo, Waa{0, k) = 0, and this proves, in particular, that 

VKaa{0,p^) = and that, in the limit |A| oo, 9^Waa(0,p^) = 0. 

Using that is invariant under the action of Ti, we see that the third of (|3.42p implies 
that (1 - e*P^(^i-^=^))Wab(A;o,P^) = 0. Since e*'^^^!-^^) = _g»^^/3 ^ ^^^^ identity proves, in 
particular, that Wa(,(0,p^) = 0, and dkoWabiO,!^) — 0. This concludes the proof of item (ii). 

Now, combining the first of (|3.43p with the first of (|3.46p . we find that Waa{ko,k) = 
Waai—ko,k)*, which impUcs, in particular, that Re{9fe(, Waa(0,p*^)} — 0. 

Finally, let Wa6(0,p^ + k') ~ aik'i + a2k'2, modulo higher order terms in k' . Using that 
^ - ^ "l/^ V3/2\ ^^.^^ ^^^^ 

v3/2 -1/2 J 



'^k[ + = e-*""'/3 \ai{k'i/2 - V3fc^/2) + a'^{V^k[/2 + fc^/2) 



(3.48) 



which implies = —iuja2- Moreover, using the third of (|3.43p we find that c 
and using the third of p.45p we find that a"^ — —ol^^ . Therefore, — 



-{a 



2 



and we see that ctj is real and odd in w, that is 03 ^ fo'' some real constant a. Therefore, 
ai = —iLoa2 = —ia, and this concludes the proof of item (iii). ■ 



C. Free energy: The infrared integration 

Multiscale analysis. In order to compute (|3.33p we shallproceed in an iterative fashion, using 
standard functional Rcnormalization Group methods [l^ [3]. As a starting point, it is 
convenient to decompose the infrared propagator as: 



("•)(x,y) = ^ e-^'^-(-^-^~)#o)(x,y) , (3.49) 
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where, if k' = (/cq, k'), 



-1 



5i^°)(x,y) = -^ E Xo(|k'|)e--'(--) ( (3.50) 

and =Vpx Vt, with - {^fh + ^62 - . [t^] + 1 < "i."2 < [|]}. 

Correspondingly, we rewrite ^(* '" ) as a sum of two independent Grassmann fields: 



= E ^''"""'"^i^i'i. (3-51) 

and we rewrite (|3.33p in the form: 

S,,^ = e-^l^l^» I P,„,^„(dvI;(<o))e-v-'(*<=^")) ^ (3.52) 

where V(°)(*(-°)) is equal to V(*(* '' )), once ) is rewritten as in ((33T|) . i.e., 

V(o)(^'(^o)) = (3.53) 

00 LUi,... ,a;27i — i n 

-y(/3|Ar2" V V V [TT*L^"^^ 

A^^"^' [11 k2j_i,0"j,P23-l,1^2j-l k2j,a-j.P2j,(^2j _ 

n=l o-i,...,CT„ = ti pi ,...,P2n = l,2 kj,...,k2„ J = l 

2n 

■W^i°)^,^(k;, . . . , k^„_i) <5( E(-1F(Pf' + k;-)) = 

00 ^ n 

n— 1 J — 1 
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with: 

1) w = {uJi, . . . ,uj2n), 21= (cri,...,cr„) and p*^ = {0,pp); 

2) W^ZJK, Kn-l) = W^nAK + Pp' : ■ ■ ■ , Ku-1 + Pf ^66 Q ; 

3) the kernels p.uji^i^ ■ ■ ■ j^2n) are defined as: 



(0) 



,X2„) 



(3.54) 



Moreover, P-^g^Aoid"^^^^^) is defined as 



,(<o)- 

.cr.p,uj 



Xo(|k'|)>0 

n n d^it'iJ^i^ 

Xo(|k'|)>0 

exp{-(/3|A|)-i Y: E Xo \|k'|)*g°Vo,.(k')^g?:j , 



(3.55) 



where: 



uj=±,<7=U k'G-Dg 



-v{k'+p^) -iko 

-iCo^o + so(k') co{ik[ - tjfc^) + to.c^(k') 

co(-^fci - wfcy + i5 ^(k') -^Cofco + so(k') 



A/q is chosen in such a way that / P^^ Aoid"^^^"^) = 1) Co = 1, cq ~ 3/2, sq = and \tQ cj(k')| < 

C|k'|2. 

It is apparent that the v]/(<o) field has zero mass (i.e., its propagator decays polynomially at 
large distances in x-space). Therefore, its integration requires an infrared multiscale analysis. 
We consider the scaling parameter 7 > 1 introduced above, see the lines preceding p.28p . 
and we define a sequence of geometrically decreasing momentum scales 7'', h — 0, ~2, . . . 
Correspondingly we introduce compact support functions //i(k') = xo(7 ''|k'|) — Xo(7 ''^^|k'|) 
and we rewrite 



Xo(|k'|)= M^') 



(3.56) 



h — ~OD 



The purpose is to perform the integration of (|3.52p in an iterative way. We step by step 
decompose the propagator into a sum of two propagators, the first supported on momenta 
~ 7'', /i < 0, the second supported on momenta smaller than 7''. Correspondingly we rewrite 
the Grassmann field as a sum of two independent fields: v]/(<'») = vj/C*) -(- \[/(<'»-i) and we 
integrate the field vj/C*). In this way we inductively prove that, for any h < 0, (|3.52p can be 
rewritten as 



S,,^ = e-^I^I^Wp,,,^Jd*(^''))e-^'^'(*'^'^') 



(3.57) 



where Fh,AhM''^ will be defined recursively, Xh{\k'\) = EL-co /fc(k') and P^.^aM'^^-'''^) 
is defined in the same way as P^g^Aoid'^^-^^) with Xo, ^o,w, Co, cq, so, io,w replaced by 

Xh, Ah,uj, Ch, Ch, Sh,th,uj, respectively. Moreover V'-'^^O) = and 



v(^)(*)=5:(/3iAi)-^"E y: [n 



(<h) + 



i<hy 



(3.58) 
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Note that the field ^'k^^'. whose propagator is given by x/i(|k'|)[A£''''(k')] ^, has the same 



support as x/i, that is on a neighborood of size 7 around the singularity k' = (that, in the 
original variables, corresponds to the Dirac point k = p"^). It is important for the following 

to think W^2,'j'pwJ h<Q, as functions of the variables {C,kiCk\h<k<o- The iterative construction 
below will inductively imply that the dependence on these variables is well defined. 
The iteration will continue up to the scale where /i^ is the largest scale such that 



007"''^' < ^Qh, (3.59) 

where ao is the constant appearing in the definition of xo(|k'|). By the properties of C,h that 
will be described and proved below, it will turn out that hfj is finite and larger than log^ 2a^' 
The result of the last iteration will be j^, i.e., the value of the partition function. 

Localization and renormalization. In order to inductively prove (|3.57[) we write 
where 

^^''' = m E E "E"*gt!rp..^g'.L,.<;!(..)(k^ (s-sd 

'"^l I <Ti,0-2=Ti P1,P2=1,2 k' 



and 7^V('*) is given by (|3.58p with Yl^=i replaced by Yl°^=2^ that is it contains only the mono- 
mials with more than four fields. 

Note that in p.6ip the w-index of the fields is the same; this follows from the fact that 
in the terms with different lo's the momenta verify k'j^ — kj + p"^ — t^f^ = "-i^i + "2^21 for 
some choice of rti, rt2, and such a condition cannot be verified if k'j^, k2 are in the support of the 
\^{<h) fields, because p"^ — p]^" ^ A* and 2ao7 is smaller than 47r/3 — 47r/(3'\/3), see the lines 
preceding p.28[) and the discussion after (|3.40p . 



Remark. The fact that the quadratic terms with different w's, i.e., the one particle umklapp 
processes, do not contribute to the infrared effective potential is a crucial fact, which reduces the 
number of relevant running coupling constants and, in particular, tells us that the interaction 
does not generate mass terms. Note, in fact, that the presence of one particle umklapp terms 
with a non zero contribution at the Fermi points could produce an exponential decay of the 
interacting correlations. 



The symmetries of the action, listed in Lemma 1, which are preserved by the iterative integra- 
on 
and 



tion procedure, imply that, in the zero temperature and thermodynamic limit, W^'^^ uj)i^) — ^ 



2,(pi,p2),(i^,w)^ ^ \Sh{~ik[ - ujk'2) -izhko ' ^ ^ 



Pl:P2 



for suitable real constants Zh,Sh- The proof of (|3.62p is completely analogous to the proof of 
Lemma 2 and will not be repeated here. 

Once that the above definitions are given, we can describe our iterative integration procedure 
for h < 0. We start from (|3.57p and we rewrite it as 

J F^^,^^(d.p(^''))e-^^"''(*'-''')-^^"''(*'-''')-^|A|^'^ , (3.63) 
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with 



Xh(|k'|)>0 

£V(''H*(^'')) = (/3|A|)-i^ ■ (3.64) 

Shi^ik'i - Lok2) + ^(k') -izhko + (T,j(k') 



Then we include CV^^^ in the fcrmionic integration, so obtaining 

' p^^_^_^^d^{<h))^-nV<-H^^^^')~m{Fn+en) ^ (3.65) 



where e/j is a constant that takes into account the change in the normahzation factor of the 
measure and 

-J (\,'\^( -</i-ifco + s/i-i(k') Ch-i{ik[-ujk'^)+th-i,u>{\<^)\ , , 

"^''^^ ' {c,-iHk[-cjk',)+ti_,^jk') -*a_ifco+s,-i(k') ; ^^-^^^ 

with: 

4_i(k') - a + ZhXhO^') , Q,_i(k') ^Ch + ShXhO^') , 

Sh-lik') = Shik') + ahik')xhik') , i,.-i,^(k') = t^^^k') + T;,,^(k')xh(k') . (3.67) 

Now we can perform the integration of the field. We rewrite the Grassmann field \|/(-'') 
as a sum of two independent Grassmann fields \I>(^''^i) + vj/C*) and correspondingly we rewrite 
as 

where 

4 tU\ - f -iCh-iko + Sh-i{k') Ch-i{ik[ - ujk'2) + th^i,^{k')\ 

^''-''^^^>-[chM-tk[-ojk'2)+tl_,Jk') -tCh-iko + Sh^k') ) ^'^■^^> 

with: 

Ch-^l = Ch + Zh , Ch-l = Ch + 6h , 

Sft_i(k') - Sh{k') + ah{k') , i,._i,^(k') = th,u{k') + r^,^(k') . (3.70) 
The single scale propagator is 

where 

5('')(xi,X2) = ^ J2 e-'^'(''--^)A(k')[A,_i,^(k')]"'. (3.72) 

After the integration of the field on scale h we are left with an integral involving the fields 
^(<'»-i) and the new effective interaction V'-'* ^\ defined as 



It is easy to see that V^'' is of the form (j3.58p and that Fh-i — Fh + + en- ft is sufficient 
to use the well known identity 

eh + V(''-i)(*(^"-i)) = l(-l)»+if,f (7ev(''^(*(^''-i) + ^C')); n) , (3.74) 

n>l 
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where n) is the truncated expectation of order n w.r.t. the propagator gif* , 

which is the analogue of (I3.34p with \|/(" '' ) replaced by ^'C'' and with replaced 

Note that the above procedure allows us to write the effective renormalizations Vh = (C,h, Ch), 
/i < 0, in terms of Vk, h < k < 0, namely Vh-i — Ph{vh, ■ • ■ , vtf), where [3h is the so-called Beta 
function. 

Tree expansion for the effective potentials. An iterative implementation of (j3.74p leads to a 
representation of V''*' (4''-''-') in terms of a tree expansion, defined as follows. 




FIG. 1: A tree r £ Th,„ with its scale labels. 



1) Let us consider the family of all trees which can be constructed by joining a point r, the 
root, with an ordered set of n > 1 points, the endpoints of the unlabeled tree, so that r is not a 
branching point, n will be called the order of the unlabeled tree and the branching points will 
be called the non trivial vertices. The unlabeled trees are partially ordered from the root to 
the endpoints in the natural way; we shall use the symbol < to denote the partial order. Two 
unlabeled trees are identified if they can be superposed by a suitable continuous deformation, 
so that the endpoints with the same index coincide. It is then easy to see that the number of 
unlabeled trees with n end-points is bounded by 4". We shall also consider the labelled trees 
(to be called simply trees in the following); they are defined by associating some labels with 
the unlabelled trees, as explained in the following items. 

2) We associate a label h < —1 with the root and we denote T^^n the corresponding set of 
labeled trees with n endpoints. Moreover, we introduce a family of vertical lines, labeled by an 
integer taking values in [h, 1], and we represent any tree r g 'Th,n so that, if v is an endpoint or 
a non trivial vertex, it is contained in a vertical line with index > h, to be called the scale 
of V, while the root r is on the line with index h. In general, the tree will intersect the vertical 
lines in set of points different from the root, the endpoints and the branching points; these 
points will be called trivial vertices. The set of the vertices will be the union of the endpoints, 
of the trivial vertices and of the non trivial vertices; note that the root is not a vertex. Every 
vertex u of a tree will be associated to its scale label hy, defined, as above, as the label of the 
vertical line whom v belongs to. Note that, if vi and V2 are two vertices and vi < V2, then 

3) There is only one vertex immediately following the root, which will be denoted vq and 
cannot be an endpoint; its scale is h + 1. 

4) Given a vertex u of r g 7^ „ that is not an endpoint, we can consider the subtrees of r with 
root V, which correspond to the connected components of the restriction of t to the vertices 
w > V. If a subtree with root v contains only v and an endpoint on scale hy + l, it will be called 
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a trivial subtree. 

5) With each endpoint v we associate one of the monomials with four or more Grassmann 
fields contributing to 7?,V^°' (vE'^-'*"""'^'), corresponding to the terms with n > 2 in the r.h.s. of 
p.53p (with $(-0) replaced by \I>(^:''"~i)) and a set Xy of space-time points (the corresponding 
integration variables in the x-space representation). 

6) We introduce a field label f to distinguish the field variables appearing in the terms associated 
with the endpoints as in item 3); the set of field labels associated with the endpoint v will 
be called 1^; note that is the order of the monomial contributing to V'-^-' (^''■-''"^^-') and 
associated to v. Analogously, if v is not an endpoint, we shall call 1^ the set of field labels 
associated with the endpoints following the vertex v; x(/), e(/), cr(/), p{f) and w(/) will 
denote the space-time point, the e index, the a index, the p index and the co index, respectively, 
of the Grassmann field variable with label /. 

In terms of these trees, the effective potential V^^\ h < —1, can be written as 

oo 

V('')(*(^''))+/3|A|efc+i = ^ V('')(r,*(^'')) , (3.75) 

where, if vq is the first vertex of t and ri, . . . , Ts (s = Sy^) are the subtrees of r with root wq, 
V'-'''(r, ^'(-''^) is defined inductively as follows: 

i) if s > 1, then 

V('')(r, ^ L^fJ;^ [v('^+i)(Ti, vI/(<"+i)); . . . ; V(''+I)(r„ *(<"+!))] , (3.76) 

where V(''+^^(ri, \I/(^^+i)) is equal to 7?.V'^''+^) (n, v]/(<''+i)) if the subtree contains more than 
one end-point, or if it contains one end-point but it is not a trivial subtree; it is equal to 
7^V(°)(ri, ^'(<''+i)) if T, is a trivial subtree; 

ii) if s = 1, then v('')(t,*^-''^) is equal to ^^f^ J7^V(''+l)(rl, *(^''+i))] if n is not a trivial 

subtree; it is equal to £^^^ [nV^°\¥-''+^^) - 7^V(°^*^-''^)] if n is a trivial subtree. 

Using its inductive definition, the right hand side of p.75p can be further expanded, and in 
order to describe the resulting expansion we need some more definitions. 

We associate with any vertex v of the tree a subset Py of /„, the external fields of v. These 
subsets must satisfy various constraints. First of all, if v is not an endpoint and fi, . . . are 
the Sy > 1 vertices immediately following it, then Py C UiPy.; if v is an endpoint, Py = ly. 
If V is not an endpoint, we shall denote by Qy. the intersection of Py and Py-; this definition 
implies that Py — UiQy^. The union ly of the subsets f^. \ Qy^ is, by definition, the set of the 
internal fields ofv, and is non empty if Si, > 1. Given r g 7^^„, there are many possible choices 
of the subsets Py, v G t, compatible with all the constraints. We shall denote Vt the family of 
all these choices and P the elements of Vt- 

With these definitions, we can rewrite V^'*-'(r, \|/(-'')) in the r.h.s. of (|3.75p as: 

V('')(r,*(^'')) = J2 ^^''\t,P) , 

V('')(r,P) I d^y„¥^'^HPyM'^'\^y,) , (3.77) 

where 

$^^'')(p„)^ n*if;t\'Ap(/)M/) (3-78) 

and K^'^p^\'x.yg) is defined inductively by the equation, valid for any v € t which is not an 
endpoint, 

<P (x.) - ^ nK"+')(x.J] SU^^'-^Hp., \ . . . , *(''")(P..„ \ Q.J] , (3.79) 

l—l 

where '^^^''\Py. \ Qvi) has a definition similar to (|3.78p . Moreover, if Vi is an endpoint 
is equal to one of the kernels of the monomials contributing to T^V'^'^^ (^'^-^^^), 
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corresponding to the terms with n > 2 in the r.h.s. of (|3.53p (with v]/(<") replaced by \I>(^''")); 
if Vi is not an endpoint, xiif"^^-' — K^'^'^,^\ where Pi = {Pw,w G Ti}. 

(|3.75p - (|3.79p is not the final form of our expansion; we further decompose V^''-'(t, P), by using 
the following representation of the truncated expectation in the r.h.s. of (|3.79p . Let us put 
s = Sj,, Pi = Pvi\Qvi] moreover we order in an arbitrary way the sets P^ = {/ G Pi, e{f) = ±}, 
we call their elements and we define x^*) = U^gp-x(/), y^*) = Uygp+x(/), x^ ~ :^{f~-), 
Yij — x(/,j^). Note that J2i=i l^i l ~ l-^i^l = otherwise the truncated expectation 

vanishes. A couple I = ifij,f^ji) = (fr ^fi~) '^i^ t)e called a line joining the fields with labels 
fijJvj'^ sector indices ^ i^{f^), = '^(/;^), p-indices = pif^), Pt = Pift)^ ^nd spin 
indices (T;~ = a{f^), af = a{f^), connecting the points x; = Xij and y; = y^'j', the endpoints 
of I. Moreover, if cj^" — cu'^ , we shall put wi = uj^^ = uj^ . Then, we use the Brydges-Battle- 
Federbush formula (e.g., see [H, [3) saying that, up to a sign, if s > 1, 

4^(*('''(Pi), . . . , ^^'^Hp.)) = E n Is^^^' - y'^lr.P^ / '^^^(^^ det G'''^(t) , 

(3.80) 

where T is a set of lines forming an anchored tree graph between the clusters of points x^*-* Uy^*'', 
that is T is a set of lines, which becomes a tree graph if one identifies all the points in the same 
cluster. Moreover t = {tn' G [0, 1], 1 < < s}, dPrit) is a probability measure with support 
on a set of t such that tn' = Ui ■ u^/ for some family of vectors G M'* of unit norm. Finally 
G'''^(t) is a (n — s + 1) X {n — s + 1) matrix, whose elements are given by 



h,T 



tii'd - +6 - + 



(3.81) 



with {fij,fl^ji) not belonging to T. In the following we shall use (|3.78p even for s — 1, 
when T is empty, by interpreting the r.h.s. as equal to 1, if |Pi| — 0, otherwise as equal to 
dctG"' =£,y($('')(Fi)). 

Remark. It is crucial to note that G'''"^ is a Gram matrix, i.e., defining e+ = = (1,0) 
and e_ = = (0, 1), the matrix elements in p.8ip can be written in terms of scalar products: 



where 



= (ui (g) e^- (g) e^- (g) A{-Xij - •) , u^- (g) (g 



— ik'x 



i?(x) 



E 



-ik'x 



(3.82) 



- •) = (fa,g/3) 



(3.83) 



The symbol (•, •) denotes the inner product, i.e., 

(ui g) g) Go- g) A(x - •), Ui' g) e^i g) e^' g) B{x' - •)) = 
= {ui ■ u,,) (e^ • e^/) (e^ • e„,) ■ J dzA*{x ~ z)B{x' - z) , 



(3.84) 



and the vectors fa,g^ with a,f3 = l,...,n — s + 1 are implicitely defined by p.82p . The 
usefulness of the representation (|3.82p is that, by the Gram-Hadamard inequality (see, e.g., 
[H), |det(f„,g;3)| < njl/o|lll5a||. In our case, ||f„|| < Cj^'^/^ and ||g„|| < C^f^/^ 
Therefore, ||/a|| \\ga\\ < C^f^^ , uniformly in a, so that the Gram determinant can be bounded 

by C'n-''+l-y2?i(rt-s+l)^ 
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If we apply the expansion (|3.80p in each vertex of t different from the endpoints, we get an 
expression of the form 



(3.85) 



where T is a special family of graphs on the set of points x^^ , obtained by putting together 
an anchored tree graph Ty for each non trivial vertex v. Note that any graph T G T becomes 
a tree graph on x^j^, if one identifies all the points in the sets x„, with v an endpoint. Given 
T G Th^n and the labels P,T, calling w*, . . . the endpoints of r and putting hi = /ii,*, the 

explicit representation of p ^ (x^,^ ) in (|3.85p is 



WV,p,t(xi,J 



(3.86) 



not e.p. 



Pi -pi 



Analyticity of the effective potentials. The tree expansion described above allows us to ex- 
press the effective potential V'-'*-' in terms of the running coupling constants Ch,Ch and of the 
renormalization functions crA;(k), tfc_i^(k). 

The next goal will be the proof of the following result. 

Theorem 2. There exists a constants J/q > such that, if \U\ < Uq, then the kernels 
Jxi,...,X2;) in i3. 58]) . h < —1, are analytic functions of U, satisfying, for any 



< 9 < 1 and a suitable constant C > 0, the following estimates: 



1 



dxi • • • dx2i|t^4';^,^(xi, . . . ,X2,)| < 7"^""''+'^ {C |(7|)™"-(i''-i) 



(3.87) 



Moreover, the constants e^ and Ch defined by iS. 68\) and | j*. 73^ are also analytic functions of 
U in the domain \U\ < Uq, and there they satisfy the estimate \eh\ + \eh\ < CUqj'^^^'^^'^ . 

Remark. The above result immediately implies the analyticity of the specific free en- 
ergy ffj{U) and of its zero temperature limit e{U), i.e., of the specific ground state energy. In 
fact, by construction, ffj{U) = Fq + Y^H^hpi^h + e/i), with Fo an analytic function of U, see 
the discussion after (|3.36p and in Appendix [BJ Therefore, Theorem 2 implies the part of the 
statement of Theorem 1 concerning the free energy and the ground state energy. For the proof 
of analyticity of the Schwinger functions, see next Section. 



Proof of Theorem 2. Let us preliminarily assume that, for h < —1, and for suitable 
constants c, c„, the corrections Zh,5h, cr/i(k') and T/i(k') defined in p.62p and p.64p . satisfy the 
following estimates: 

max{|^,,|,|5;,|} <c|C/|7^\ (3.88) 
sup {\\dl,aumi\\dl>r,,Am} < c„\Uh<^'+'>-^'^ . 

|k'|s.t.Xh(k')5^0 

Using (j3.88p we inductively see that the running coupling functions (^/j, c/i, s/i(k') and t/i(k') 
satisfy similar estimates: 

max {\Ch - 1|, \ch - 3/2|} < c\V\ , (3.89) 
sup {\\dl,s^{\l)\\,\\dl.t^^M\\} < c„|^|7''+'-"^' • 

\k'\s.t.Xh{k')^0 

Now, using the definition of g^\ see p.72p and (|3.66p . we get, after integration by parts, for 
any > 0, 

(2+7i}h 

II [dX9:^!9i'H^l,^2)]^J < , (3.90) 
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where n = no + ni+n2 >0 and Cn.ti is a suitable constant. 

Using the tree expansion described above and, in particular, Eqs. p.TSp . (|3.77p . p.85p and 
(|3.86p . we find that the l.h.s. of (j3.87|) can be bounded from above by 



E E E E/ n<--y' 



i=l 



(3.91) 



n ^max|detG''-^"(t.)| J] | |5i':Hx; - y,) | 



_ V not c.p. 



where 1 1 • 1 1 is the spectral norm and where T* is a tree graph obtained from T = UyTy, by adding 
in a suitable (obvious) way, for each endpoint v* , i — 1, . . . ,n, one or more lines connecting the 
space-time points belonging to x^,. . 

A standard application of Gram~Hadamard inequality, combined with the dimensional bound 
on ^^''''(x) given by p.90p . see the remark after p.8ip . implies that 

idetG'-'^^'(t.)i < cI:::j^".i-i^^'|-2(^"-i) . 1-2(^-1)) . (3.92) 

By the decay properties of g^\-x.) given by (|3.90p . it also follows that 



n ^/n^(^'-y')iifff;;^^N-yoii<c" n ^t-^-^^--^^ (3.93) 

V not c.p. ^ ^^T^, V not c.p. ^ 

The bound (|3.37p on the kernels produced by the ultraviolet integration implies that 

/n n 
i^m^\ . , m ^ — 1 ; — 1 



;eT*\u„T„ 



1=1 



where pi = l^?;*!- Combining the previous bounds, we find that p.9ip can be bounded above 

by 

EE E E^l n ^7'^"(^-l"-l-l""l-^^--^^)][n^^'''|C^|--^] (3.95) 



n>lTeTu„ P£Tt TgT u not c.p. 

- • |P,,„I=2; 



Let us define n{v) = 'Ylii-v*>v ^ ^ ^^e number of endpoints following on t and v' as the vertex 
immediately preceding v on r. Recalling that is the number of field labels associated to the 
endpoints following v on r (note that |/^| > An[v)) and using that 



E [(Ei^^^.i)-i^^ 



V not e.p. i—\ 



l-^'uo I l-ffo I J 



^ (s„-l) = n-l, 

V not c.p. 

Si, 

^ (/i.-M[(Ei^''.i)-i^''i] = E (/^.-/^.od^i-i^.i), 

V not c.p. z— 1 u not c.p. 

(ft.„ - /i)(si, - 1) = ^ [K - h^,){n{v) - I) , 

V not c.p. t; not c.p. 

we find that p.95p can be bounded above by 

EE EE cv^^"'^""'^'^""'"^"^ • 

' |P«ol=2! 

[ n i 



(3.96) 



n>lTeTh„ PS^T TeT 



iK-h^,)(3-\P^\ 



V not c.p. 



i=l 
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Finally, let n{v) be the number of endpoints following v but not following any of the vertices 
w > V and let p(v) be the number of field labels associated to endpoints following v but not 
following any of the vertices w > v. Using the identities 

V not c.p. V not c.p. 

^M/.ol Yl ^ih.-h^')\i.\ = Y[ l^^^^^'K (3.98) 

V not c.p. V not c.p. 

we obtain 



^•(3-|P.ol) 



-i- /dxi...dx2d<i,jxi,...,x,oi<E E E E^V'^ 

1 " 



V not c.p. not c.p. 



Note that, if v is not an endpoint, |Pu | — 3 > 1 by the definition of TZ. Moreover p{v) — 3n{v) > 
and jjQj g p (j){v) — 3n{v)) > n; in particular, this means that there exists at least one vertex 
V* that is not an endpoint, such that p{v*) — Sn{v*) > 1. Therefore, we get 

V not c.p. 

with the highest scale label of the tree. Now, note that the number of terms in X^tgt 

be bounded by C" Ht, not c.p. ^vl- Using also that - 3 > 1 and - 3 > |Fi,|/4, we find 

that the l.h.s. of (|3.99p can be bounded as 

i^|dxi...dx2d<:yxi,...,X20|<7''^'-|''-I^E^" E V'-- (3.101) 



/3|A| 



n>l reTh^n 



not c.p. Pev-r- V not c.p. i=l 

Now, the sum over P can be bounded using the following combinatorial inequality (see for 
instance §A6.1 of [i3|): let {pv,v G t}, with r £ Th^n, a set of integers such that p„ < J2i=i Pvi 
for all w g T which are not endpoints; then, if a > 0, 

n E7"""^^c 

V not e.p. Pt, 



This implies that 



pg'Pt V not c.p. 



Finally, using that 7''* n„ „ot e.p. 7^^^''"^''"'' < T*"*' and that, for < 61 < 1, 

TGTh.n i" not c.p. 

as it follows by the fact that the number of non trivial vertices in r is smaller than n — 1 and 
that the number of trees in 7^,„ is bounded by const", and collecting all the previous bounds, 
we obtain 



n>l 
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which is the desired result. 

It remains to prove the assumption p.88p . We proceed by induction. The assumption is 
vahd for ft, = 0, as it follows by (|3.37p and by the discussion in Appendix |B] Now, assume that 
p.88p is valid for all h > k + 1, and let us prove it for k — 1. The functions —izkko + o'fe(k') 

and 6k{ik'i — cjfcy + Tk.ui^) admit a representation in terms of w!^^^ y)- particular, 

max{|zfc|,|4|} '^xi'^X2|x-y||1^2('2(^^^)(x,y)| , (3.103) 

and 

sup {\\dl,au{y)Mdl,Tk,^m\} < Cj"^-^ [ dxidx2|x-y|"+2|T4^W (x,y)| . 

(3.104) 

The same proof leading to p.l02p shows that the r.h.s. of p.l03p can be bounded by the r.h.s. 
of (|3.102p times 7"''' (that is the dimensional estimate for |x — y|), and that the r.h.s. of (|3.103p 
can be bounded by the r.h.s. of (|3.102p times -^2/c^-(ri+2)/c (^^^j^gj-g ^-fe("+2) jg ^]^q dimensional 
estimate for |x — y|"^^). This concludes the proof of Theorem 2. ■ 



D. The two point Schwinger function 

In this section we describe how to modify the expansion for the free energy described in 
previous sections in order to compute the Schwinger functions at distinct space-time points. 
For simplicity, we shall restrict our attention to the case of the two point Schwinger function. 
The general case can be worked out along the same lines. 

The Schwinger functions can be derived from the generating function defined as 

>V((^) = log / P(d*)e-^('^)+/ '*-['^^,-,P*>^,-.P+*i,-.P^^,-.p] (3.105) 



where summation over repeated indices is understood and the variables 0x a- p Grassmann 
variables, anticommuting among themselves and with the variables Vt^ o- p. The two-point 

def 

Schwinger function S'(x — y)p,p' = >5'2(x, cr, — , p; y, cr, +, p') is given by 

5(x - y)p,p^ = ^1- _ ■ (3.106) 



0=0 



We start by studying the generating function and, in analogy with the procedure described 
before, we begin by decomposing the field in an ultraviolet and an infrared component: 
^ ^ -qjin.v.) _^ ^(t.r.)^ ^-^^Yi = J2uj=± e'^^^'x-i'^p'lL. After the integration of the " ) 

variables, and after rewriting 0x,cr,p — Sw=± ^^^^'f'x.cr,p^ S^t: 

gW(0) ^ g-/3|A|Fo+5(^°)(0) J p^^^^^(dvI/(<0)) . (3.107) 

where 5*^-"^ ((/>) (chosen in such a way that 5''--°-' (0) = 0) collects the terms depending on 4> but 
not on \l>(-0) and B^"^ (^''-''^ </)) the terms depending both on and ^'(-O) generated by the 
ultraviolet integration. 

Proceeding as in Section IIII C[ we inductively show (see below for details) that, if /i < 0, 
gW(</)) gg^j-^ ^g i-g-vifritten as: 
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where: J dk' must be interpreted as equal to (/3|A|) ^ X^keP" ! ij('0(\i>(<'»)^ 0) can be written 
as + with W^"^ containing the terms of third or higher order in (j) and 

bII'\^^^^^) of the form 



^ ^ ■ .(T.po.UJ ■ .(T 



+'?'-,<Ti,pi,a;i * (^a;i,pi,p2 * „ (</i)+ ;^^(</i)- J * '-^a;2 ,P3 ,P4 * ,<T2 ,P4 ,^2 

with 

GSptpkk')^ ^ 5Sip"(k')QLt,P",P' (3-110) 
and p p' defined inductively by the relations 

*^k',i^, p,p' — "^k'^cj^p^p' '^'^2,p,p",(w,w)'>*^ '''-^i^,p",p'V*' -I ' ^k',i^,p,p' — "P.P ' 

where VFj^l^.^ is the kernel of rV^*), as defined in ([SIT]). In p.llOp . gi^^ is defined as 

5«(k0=5<"--)(k' + p5^)[]l(||k'||<||k' + p^^-p^-||) + i]l(||k'|| = ||k' + p^-p^"| 

where p"^ '''^ (0,p^). Note that, by the compact support properties of ^i'*^(k'), if ^^''^(k') ^ 0, 
/i< 0, then .g(^'(k) = for \ j -h\>l, so that 

'^k',CJ,p,p' — ^^2,p.,pi,(u,i^)\'^ )yuJ.,puP2\'^ )^k\L0,p2,p' ' 

and, therefore, proceeding by induction, we see that on the support of gi''''(k') we have 

lid - 111 < C^IC^It'" , W^'QiVl < C„|[/|7(^-")'' . (3.112) 

In order to derive (|3.112p . we used Theorem 2 and the decay bounds (|3.90p . 
Using (|3.112|) . the definition (|3.110p and the decay bounds ((3J0)) . we find that 

dx|xn|G('')(x)||<Q7-(^+^')^ (3.113) 

Let us now prove p.lOSp . We proceed by induction. For h = p.lOSp is clearly true (it 
coincides with (|3.107p ). Assuming inductively that the representation (|3.108p is valid up to 
a certain value of h < 0, we can show that the same representation is valid for h — 1. In 
fact, we can rewrite the term V^*) in the exponent of (|3.108p as V^''^ = /IV^''^ + TZV^''\ as in 
p.60p . and we can "absorb" CV^'^'^ in the fermionic integration, as explained in Section fill C[ 
see p.63p - p.65p . Similarly we rewrite 

- V^'^U'if^-'''')^ f dyW^''^ w ^^x v)*^-''^T + - TZV^'^H^^-'''') (3 114) 

a,T,(<'')± ^ ' ^■>'^*^2.Xp..p'),{uj,uj)\^'y)^y,a,p',uj^ i<h)±"^'^ {O.ll'i) 

(y^X,Cr,p,LJ •> Cyr'x,Cr,p,LJ 



This rewriting induces a decomposition of the first line of p.l09p into two pieces, the first 
proportional to w'^^ ^ the second identical to the first line of p.l09p itself, with V''*' replaced 
by 7^V(''^ that we wiU caU 7^B^''' (*(-'')). We choose to "absorb" the term proportional to 
w'^'' into the definition of Q^'*', and this gives the recursion relation (|3.11ip . Moreover, note 
that combining 7^B^''^(*(^'')) with 7^V('')(*(^'')) wc find: 

7^v('')(**-''^) + 7^B^''^(*(^'')) = 7^v('')(**-''^ + * </>) + w^^l , (3.115) 
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with wll^l containing terms of third or higher order in (j). We define W^"* — W^^'' + W^^l . 
After these splittings and redefinitions, we can rewrite (|3.f OSp as 

^WW ^ ^-P\A\iF,+e,) + S'>--H'}>) J p^^_^ ,i,^_^(dM,(<'-l)) J (3.116) 

Integrating the field vJ/C'), we get the analogue of p.73p : 

^ g^e,/3|A|-V('^-i)(*(2.-i)+G(^)*^)+/dk'^+Q<';)g<'')(kOQ<';'0^,-<':-^ . (g^^^^) 

with G'-'*-' defined by the recursion relation (|3.110p and ^ term of third or higher order 

in <j). Eq. (j3.117|) can be proved by making use of a formal change of Grassmann variables 
ii'^, = ^f^, - g('*Hk')QL'?^0k', as described in Ch.4 of 0. At this point it is straightforward to 
check that the final expression for e^'^'^^ that we end up with is given by the r.h.s. of p.lOSp . 
with h replaced by /i — 1, and the inductive assumption is proved. 

From the definitions and the construction above, we get 

5p,p.(x - y) = ^ e-*J^(-~-«~)5^,p,p,(x -y)^Y. ^-'^^^^-^ ■ (3.118) 

o;— ± a; — ± 

h— — OQ 

Analyticity of 5p_p'(x — y) follows from this representation and the results of Theorem 2. 
Concerning the representation p.lOp . let us take the Fourier transform of Su,p,p'{x — y). If we 

define /ik = min{/i : gi''^(k') ^ 0}, we get, for k' inside the support of ^[^^J p ^ 

hk + l 

Juj,p.,p'\i^ ) — Vk',w,p,piyt^,pi,P2v"- ^k',w,p2.p' 
j=iik 

- E G^^LA^')<XpZ.U^.^)^^^^^^ ' (3.119) 

j=hk 

which readily implies (j2.10p : in fact, using the explicit expression of gi''^ and the inductive 
bounds on see p.ll2p . it is easy to see that the term in the first line of p.llOp can be 

written as in (|2.10p and that their only singularity is located at k' = 0. 

The contributions from the second line can be bounded using the bounds on 14^2''^ proved 
in Theorem 2, and we find that they can be bounded by C\U\j'^'^'^^^^^\ which means that 
they only contribute to the error term appearing in (|2.10p . This also implies that no other 
singularity, besides the one at the Fermi points, can be produced by such terms. 

Finally, if k does not belong to the support of write 

5p,p'(k) = S^;:;-\k) = g("-"-)(x - y) - * I^2,(p,,p.) * 9ty)i^ ' (3-120) 

with W2,p defined by p.35p . The bounds discussed in Section lTlI Bl and Appendix [B] imply that 

S^^p]''\'x — y) decays faster than any power, so that no singularity can appear in its Fourier 
transform. 

A similar expansion can be obtained for higher order Schwinger functions, but we will not 
belabor the details here. This concludes the proof of Theorem 1. ■ 
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APPENDIX A: THE NON-INTERACTING THEORY 



In this Appendix we give some details about the computation of the Schwinger functions of 
the non interacting theory, i.e., of model (12. ip with U ~ 0. In this case the Hamiltonian of 
interest reduces to 



i=l,2,3 

with A, ai , - defined as in items (l)-(4) after (l2l1) . 

First of all, let us remind that, being i?o,A quadratic, the 2n-point Schwinger functions satisfy 
the Wick rule, i.e., 

(T{*.,,.„p, • • • *x„,.„.p, ,.;.p; • • • ^L<.P'J),^^ = - det G , 

Moreover, every n-point Schwinger function ^^'^(xi, ei, cti, pi; . . . ; x„, e„, (t„, p„) with 
Sr=i £i 7^ is identically zero. Therefore, in order to construct the whole set of Schwinger func- 
tions of Ho,\, it is enough to compute the 2-point function '^(x— y) = (T{^'~ CT.p^y o- p'})^ 
and in order to do this, it is convenient to first diagonalize -ffo.A- Let us proceed as follows. 
We identify A with the set of vectors in a fundamental cell, and we write 

A = {mai + Ti2a2 : < rii, ?^2 < i - 1} , ©3) 

with fli = i(3, VS) and 02 = ^(3,— \/3). The reciprocal lattice A* is the set of vectors such 
that e*^^ = 1, if 5" G A. A basis 61, 62 for A* can be obtained by the inversion formula: 

&2I 622 ) ^" \ ai2 0,22 



bn 612 \ _2^r" "21 I ^ 



which gives 



6i = f(l,V3), 62 = y(l,-V3). (jAl5) 



We call I?L the set of quasi-momenta k of the form 

A; = — + — ^62 , 7711,7712 eZ, (06) 

identified modulo A* ; this means that V]^ can be identified with the vectors k of the form (|2.2p 
and restricted to the first Brillouin zone: 

■Dl = {k = + : 0<77il,7772 <i-l}. (07) 

Given a periodic function / : A ^ R, its Fourier transform is defined as 
which can be inverted into 

/(fc) = E e-^^^'fix) , keVL, dug) 

where we used the identity 

E e''^^ = |A|<5,^,g (El 10) 
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and S is the periodic Kronecker delta function over A* . 

We now associate to the set of creation/annihilation operators ai , 6* - the corresponding 

set of operators in momentum space: 



Note that, using ((IT8l) - (|XrT0)) . we find that 



are fermionic creation/annihilation operators satisfying 

{«L'«1.'} = \MSj:rk'Se,-e'Sa,.' , ^ |A|%.£,<5,._.5.... (jAll3) 

and |a- ,51 | ~ 0. With these definitions, we can rewrite 

^ k.a k' ,(7' 

^o.A = - E E (4A^+.v. + - ®i4) 

u p /.^ /.^ /L^ \ k,<j k',cr k',a k.a' 

= -TTT WrO^ +vj;bi aZ ) , 

' ' keVL '^=Ti 

with 

v-^^Y. e'^'-"'^'' = 1 + Se-^i'^i cos . (015) 

i=l 

The Hamihonian -ffo.A can be diagonahzed by introducing the fermionic operators 
in terms of which we can re-write 

^o,A = ^ E E ( - \^kH.^k,^ + J ' ©17) 

with 



= y'(l + 2cos(3fci/2)cos(\/3fc2/2))V4sin^(3fci/2)cos2(V3fc2/2) , (jUlS) 
which is vanishing iff fc = p^, = ±, with 

Now. fo, S 6 A, we define of , = |A|-' and = |Ar' Ees„, 

moreover, if x = (a^o,^) we define a^.cr = e^''-'^'^°a^ ^e~^°''^'^° and /3x,<t = e-^"''^""'^±^g--ffo.Aa:o_ 
A straightforward computation, see, e.g., Appendix 1 of jll, shows that, if — /3 < xq — yo < 

(T{a_a+^,})^_^= (IS120) 



1^1 



J2 e-^'^^'-y^ l{xo - yo > 0) ^. . -t{xo-yo<0 
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(T{/3-./3+.,}),^ = (1321) 



|A| / 



{xo-Vo)\vj:\ p-{xo-yo+P)\vg\ 



Y: ^"'^'-^ [Hxo -yo> 0)^^^3^ - l{xo - 2/0 < 0) — 



fee Pi 



and (T{a,^,/3+^,})^^^ = (T{/3,_,a+^,})^^ = 0. A priori Eq.dMH) and dSEl arc defined 
only for — /3 < xq — yo < P, but we can extend them periodically over the whole real axis; the 
periodic extension of the propagator is continuous in the time variable for xq — ijq ^ and it 
has jump discontinuities at the points xq ~ yo /3Z. Note that at xq — yo = /3n, the difference 
between the right and left limits is equal to (— l)"^^ ^;, so that the propagator is discontinuous 
only at X — y = /3Z x 0. For x — y ^ /3Z x 0, we can write 

(T(o„o+,. = J„ |i j: ^ . fan, 

Note indeed that for a^o ~ J/o ^ P'^ the sums over in (jAI22p are convergent, uniformly in M; 
if xq — yo = (3n and x ^ y, the r.h.s. of (IAI22P is equal to 



lim (T{a,,,a+^,}) ^ + lim (T{a,.,a+,,}) J = (024) 



A similar remark is valid for {T{l3^ ^(3^ ^,}) ^ ^. If we now re-express ^ and (3^ ^ in terms of 
o- and o-' using (|All6p . we get (|2.6p . Note finally that if x = y 

so that the diagonal part is vanishing; on the contrary, using (|AI16|) and the fact that 

1 ^ / e^l^^ e-^l^gl _ 1 ^ 
- -2|A| ^ Vi + 6/3^1 ^ l + e-z^l^el^ " 2' "^^^^^ 

and this explains why there are no quadratic terms in V{'^), see p.l9p . 



APPENDIX B: THE ULTRAVIOLET INTEGRATION 



In order to prove Eq. p.35p - ()3.37p . a simple application of p.80p and determinant bounds is 
not enough, because (7^" '"'^(x) does not admit a Gram representation, which is a key property 
needed for the implementation of standard fermionic cluster expansion methods. As mentioned 
in Section IIII Bl a way out of this problem is to decompose the ultraviolet propagator into a 
sum of propagators, each admitting a Gram representation, and performing a simple multiscale 
analysis of the ultraviolet problem, in analogy with the standard strategy for ultraviolet prob- 
lems in fermionic Quantum Field Theories fl5| . This multiscale analysis is very similar to 
(but much simpler than) the one describe in Section IIII CI it has been performed in several 
previous papers [3, 0, and it is reported here for completeness. 

Let M be the integer introduced after l|2.6p . and let us write 

5("-''-)(x)=ygW(x), 01) 

h=l 
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where 



^^^''^W = 7¥aJ E /"-(k)i^/.(fco)e-*''5k , 02) 



fi.L 



with 7Ji(fco) = Xo(|fco|), Hhiko) = xo(7 ''^^|fco|)-Xo(7 ''^^|fco|) and Hm is the smallest integer 
such that fu,v.O^)Hj{ko) = for all j > Hm (note that Hm ^ log(Af //?)). Note that g^^'Ho) = 
and, for any integer K >0, ^(''^(x) satisfies the bound 



1 + (7'*|a;o|/3 + fIa)" 



where | • |/3 is the distance from the origin on the one dimensional torus of size [3, while | • |a is 
the distance on A. Moreover, g'^^^yC) admits a Gram representation: ^(''^(x — y) = /dz Aj|j(x — 
z) • Bft(y - z), with 

Bhi^) = ^ E V/«...(k)i/.(fco)e-*-f % ©4) 

and 

= j dz|A,.(z)|2 < ^7-"" , \\Bhf < Cj'" , 05) 

for a suitable constant C. 
Our goal is to compute 

g-/3|A|fo-v(*(")) ^ lij^ /■ p(^vl/[i^"Ml)e^(*'"^-'+*'''"*'') , P6) 

Af^oo J 

where P((i5't^''''^l) is the fermionic "Gaussian integration" associated with the propagator 

J2h=i g^'^H'^) (i-e-, it is the same as F(d*("-" ))). We perform the i ntegra tion of (|BT6)) in 
an iterative fashion, analogous to the procedure described in Section IIII CI for the infrared 
integration. We can inductively prove the analogue of p.57p . i.e., 

g-/3|A|i=^o-V(*(' '-') ^ g-/3|A|i=^;, J p(-^^[l,/i])gV('''(*(--> + *li-''l) 

where P((i4'[^'''l) is the fermionic "Gaussian integration" associated with the propagator 
Eti5^'Hk) and 

oo „ n 

v('«)(v,[i,'«i) ^Y.T. rfx,...dx,4n*ic\tiu,p.-.^i^"-^^^^ ■ 

n— 1 p,CT j — 1 

08) 

In order to inductively prove (|BI7P - (jBI8P we simply use the addition principle to rewrite 

y" P(d*[i'''l)e^"'' y J p(d5»('0)e'^<'''(*<' '' '+*'''''"''+*<''') ^ 

09) 

where P((i5'^'')) is the fermionic Gaussian integration with propagator ^'^''^(k). After the inte- 
gration of \I>('')) we define 



e 



V(h-i)(^(....)+^[i,'.-il)_^|A|^^ = / p(^^(/i))gV(''>(*(-'-'+*li''-il+*<''') ^ ^^Q^ 
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which proves (|BI7|) . In analogy with (|3.74p we have 

n>l ^' 

As described in Section IIII C[ the iterative action of can be conveniently represented in 
terms of trees r g TM-h,n, where TM;h,n is a set of labelled trees, completely analogous to the 
set Th^n described before Eq. (|3.75p . unless for the following modifications: 

1. a tree r S TM;h,n has vertices v associated with scale labels h + 1 < < /im + 1, while 
the root r has scale h; 

2. with each end-point v we associate 

y(^[i,'«iv/])^ with V{-^) defined in H^J^ . 

In terms of these trees, the effective potential V^^^ < h < Hm (with V^°\-^^' ■''■'>) identified 
with V(5''* '" -')), can be written as 

oo 

n=l reTM.h.n 

where, if wq is the first vertex of r and ti, . . . , {s — Suq) are the subtrees of r with root Wq, 
V^'^\t, vI/[1''']) is defined inductively as follows: 

i) if s > 1, then 

vW(r, vi/H.'^l) = [-pc^+i^n, . . . ; v(''+^\t,M' , (013) 

where v(''+^' (n, is equal to v(''+^' (r^, if the subtree n contains more 

than one end-point, or if it contains one end-point but it is not a trivial subtree; it is 
equal to if is a trivial subtree; 

ii) if s 1, then V^^'\tM-^^) is equal to [vC'+i) (n, if n is not a trivial 
subtree; it is equal to El^^[V{^^'^^^+^^) - ^(^'[^'''1)] if n is a trivial subtree. 

Note that, with V{'^) defined as in (|3.19|) and with the present choice of the ultraviolet cutoff 
(such that ^(''HO) =0), we get £l^^[V{^^^''^+^'^) -V{^^^'^'^)\ = 0. This implies that, if w is not 
an endpoint and n{v) is the number of endpoints following v on r, and if r has a vertex v with 
n{v) ~ 1, then its value vanishes: therefore, in the sum over the trees, we can freely impose the 
constraint that n{v) > 1 for all vertices w G r. From now on we shall assume that the trees in 
TM;h,n satisfy this constraint. 

Repeating step by step the discussion leading to (|3.77p . (|3.85p and (|3.86p . and using analogous 
definitions, we find that 

V(")(r,P) = ^ / dx„„$(^'')(P„jW^(;'j^,^(x.J ^ V('')(r,P,T) , JHM) 

where 

feP^ 

and 

not e.p. 

(11116) 

Moreover, G''"'^^(t^) is a matrix, analogous to (|3.8ip . with ^- replaced by 1 and 
replaced by g^'^\ 



Pi 'P 



JeT„ 
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h h-^1 




hw-l hM h^+l 



FIG. 2: A tree r G TM:h,n with its scale labels. 
As in the proof of Theorem 2, we get the bound 

-1- /dx,...cix,,i<i(xi,...,x,,)i<^iL/r E E E/n^( 

^' ' n>l reTM:/! „ pet't TgT ZeT* 



xi - yi) ■ 



n ^max|detG"-^"(t,)| J] | l^^^^lx; - y/)| 



not c.p. 



(Ell?) 



and, using the analogues of the estimates (|3.92p . (|3.93|1 and p.94p . taking into account the new 
scaling of the propagator, we find that (jBllZP can be bounded above by 



E E E E^"ic/r[ n ^7-'-'--) 

. . _ _„ . Oil. 



(11118) 



V not c.p. 



Using (|3.96p we find that the latter expression can be rewritten as 



E E E E ^"it^r7-"^"-1 n 



(11119) 



n>lTeTM-.h ^^^^ TeT 

|Punl = 2i 



V not c.p. 



where we remind the reader that n{v) > 1 for any t G TM-h,n- Performing the sums over T, P 
and T as in the proof of Theorem 2, we finally find 



m 

which is a special case of p.37p . The proof of the general case is completely analogous. 



(11120) 



APPENDIX C: GRAPHENE AS ASYMPTOTIC INFRARED MASSIVE QED2+1 

In this Appendix we describe the relation between 2D graphene and a regularized version of 
euclidean QED2+1 with a massive photon, massless dirac fermions and an ultraviolet cut-off. 
Let us first introduce the model of regularized QED2+1 and let us next describe its connections 
with the graphene model described in this paper. 
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We consider the following generating function for euclidean QED2+1: 

^WL.a{J,4') — / p(^dip)P{dA) e-f ''''('^o-4f,,x5/'x7M'/'x+Jf,.x5/'x7M'/'x+0xV'x+<?xi/'x) Jell) 



where: 

1. if c is the speed of light, / dx is a shorthand for a'^c"^ SxeAa' " lattice spacing 
and Aa is a periodic lattice of side Lc~^ in the time direction, of side L in the two 
spatial directions, and with sites labelled by xq = noac~^, x = na, with La~^ integer and 
Uf, = 0, . . . ,La^i - 1, = 0,1,2; 

2. summation over repeated indices = 0, 1, 2 is understood; 

3. eg is a constant, and (j) are the external fields, and 7^ are euclidean gamma matrices, 
satisfying {7^,71/} = —2(5^^, and defined as 

/Oao\ /0cr2\ /OaA 

70 = , 71 - , 72 = , 02) 

\<Jo 0/ V-CT2 0/ V-fTi 0/ 

with cr^, /X = 0, 1, 2, the Pauli matrices: 

''-(ID- ''-{ID- ---{°/o)-- ™ 

4. is a 4-components Grassmann spinor of components V'x,i, * = lj---j4; moreover, 
V'x = V'x7o, with "0+ a Grassmann spinor of components V'x ii 

5. let I? be the set of space-time momenta k with fco — 2TTcL~^{mo + i), k — 27r_L~^m, 

with = 0, 1, . . . , La~^ — 1, /i = 0, 1, 2; if we define ip^- = a'^c"^ ExeA e^*''''V'x,i, the 
fermionic integration can be written as 

1 4 

^(#) = n n^<»^^k,^] expj - ^ Xo"'(|k|)^k*K0k} , 04) 

k6'Di=l keXi 

where ^ = 7^fc^, c is the speed of light, Z is the wave function renormalization. A/" is a 
normalization constant and xo is the cut-off function introduced in Scc lIII Bl 

6- ^/j,x is a euclidean gaussian boson field associated to the gaussian measure P{dA) with 
covariance 

.,,.(x - y) = S,M^ - y) ^ 5' Yl ^-''"^""'^^^ ' 05) 

with M > 1 the "photon mass" . 
Integrating out the gaussian boson field, we can rewrite: 

gWi,a(J,0) _ J p(^(l^^Q~^'^^^+ J dxjf,_^tl)^f^ip^+ J dx{<p^^^+<l,^ip^) 06) 

where 

V(V) = " f / ^^^y (V;x7p^x)«(x - y){i^y7^^y) ■ 07) 

The four dimensional version of the above model was studied in fl^ by RG methods; the analy- 
sis (that can be repeated for the three dimensional model considered here without any relevant 
difference) is essentially identical to the one described in this paper for the 2D Hubbard model. 
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Note in particular that, identifying the spinor ip]^ with (4'k,(T,i.+, ^k,cr,2,+, ^k.(T.2.-, ^k.o-,!,-), 
both the fermionic integration P{d^p) and the effective interaction V{'ip) are invariant under a 
number of symmetries, analogous to (4)-(8) of Lemma 1, i.e., 

(4 ) ^(^/^^^ /j.^ ) ^ ^ {ko .k2 ^ki) ,(T^2,uj^ (kQ,k2,ki)^fT,2^uj ^ (ka ,ki ,k2) .(7,1^^^ 

(5') 'I'ko-pcj ~^ ^-\^ a p -uj^ ^ ~^ '^*' where c is a generic constant appearing in Pi^d^) and/or 
in V(V'); " 

(6 .a) **(fco,-fci,fc2),<T,2,w' 
(6'.b) J. J. , -> I 1 ^ ; 

[kQ,ki,k2),(7.p,uj {kQ,ki. — k2).(T,p,~i^^ 

(7') i/^ ~ iif^ 

{ko ,k) ,cr.p.Lj {kQ^~k).a,p,—uj^ 
(8') - 

^ ' (ko,k),a,p ^ ' {-ko,k),cr,p' 

It is important to note that, in addition to the symmetries (4')~(8') above, QED2+1 
also admits extra symmetries, related to its relativistic invariance, which have no counterpart 
in the Hubbard model, e.g., 

(9') lAk e^I'^'O'^il-i/j^-ij., -01^ '4'jj-i-^e~^^''°'''^^\ where Rgk — {koCosO ~ cki sm9,ki cos9 + 

c^^fco sin0, fc2)- Note that in the limit L,a^^ — > 00, there is no constraint on the choice of 6, 
while for finite L and a we are forced to choose 9 — 7r/2. The proof of the invariance of the 
model under the symmetry (9') is a simple consequence of the remark that 

g-|[7o:7i](-^ij^^^^^2)e^f'^°^'^il = (70 COS 6* - 71 sin 6' ,71 cos6' + 70 sin 6* ,72) , jOS) 

which implies that X^keP K^k is invariant under (9'). In particular, if = 7r/2, in terms of 
the components ^'k.p.o-.c^ of the spinor, (9') reads as follows: 

^(A:o,fci,fc2),o-,p,w '^^^^ + ^'^2)p,p'^{cki,-c-^ko,k2),cr,p',uj , 

%oM.k2)..<7,P.^ ^ ■^%k^,-c-^ko.k2),'ryJ'^^ + ^'^2)p',p ■ 09) 

This symmetry also implies that the kernels of the quadratic part of the effective potentials 
have a special structure. In fact, repeating the proof of Lemma 2, using symmetries (4')-(8'), 
and if W^2'(pj p^) ^0^) is the kernel of the quadratic part of the effective action at scale h, we 
find the analogue of p.62p : 

k'9k'M^i'; , (0) = f w u ^ ~!'^'^) ■ mm 

2,(pi,p2),c^v ; ySh{-iki ~ Luk2) -izhko J ^^^^ ^ 
On the other hand, for QED2+1 we also know that 

2^'^K'7.,^\5,,{-ik,-u;k2) -iZHko j^k,...,.. lyil) 



ke-D 



must be invariant under (jCI9p . which implies cz^ — i.e., the speed of light is not renormalized. 

The same proof shows that if, in relativistic notation, ipkk.^W^ipk is invariant under (4')- 
(9'), then = Cj^, for some constant C. This is precisely the same as in four dimensional 
euclidean QED. Therefore, we can repeat step by step the construction in and, in particular, 
we find that the following Ward Identity (WI) is valid: 

iZeop^ (j^^p; VkV'k-p) = e[(i/'k-pV'k-p) - (V'kV'k)] (1 + -ffo(k,p)) , 012) 

where: 

1. OM,p;^k^k-p) = /dx/ dze-»P(--y)e''^(— y)(j^,,;7/-x^y), with 

OM,z;V'xV'y = , lim f., , , - ; (E113 
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similarly, 

L'^ ^ ' L,a-i^oo OCPxOCPy J=0=0=O 

2. e = Co — c+Gq + ©(cq), with c+ a suitable constant; 

3. the correction i/o(k, p) is such that, for momenta k, p, k — p all on the same scale h (i.e, 
all belonging to the support of fh, for some finite h <0) 

\Ho{k,p)\<C\e\j"^ , 015) 

for some < 6* < 1. 

Note that the above WI differs from the formal WI obtained by neglecting the ultraviolet 
cut-off, because of the presence of the renormalized charge e = bq — c+Gq + O(eQ) and of the 
correction Hq (k, p) . 

There is a strong connection between the above model and the Hubbard model. Indeed from 
p.l20p we know that 

5(x-y) = 5(i)(x-y)+ E e'^^^^ (^~-^~)^(^°)(x - y) , (016) 

where S^^\:ii — y) is given by the sum in the second line of (|3.118p restricted to ft. < 0, and 
|S''^^^(x — y)| < C|x — y|^^^^ for |x — y| > 1; this means that, for large distances, S'^^^ is 
asymptotically negligible with respect to s!;^^\x — y). 

By p.lOp and the construction in Sections IIII CI and IIII D| we expect that the Grassmann 
spinor {'^'^^\ ^, +, ^k^^ 2 ^k^ir i pl^Y^ the same role as the spinor ^Ak in the QED2+1 

model. In order to make this intuition precise, it is convenient to combine S^^\x. — y) in the 
following matrix 



where S^^^'^ is the transpose of S'ir"''. tj(x— y) will play the same role as the correlation (V'xV'y) 
defined in (|CI15|) . in a sense to be made precise below. Similarly, the role of (j^^.z; V'xV'y) will 
be played by the correlation S2,i-p,{z; x, y), /i = 0, 1, 2, defined as 

52,i;M(z;x,y)p,p, = {T{n^„,p^l,y^+„,.<J,,^,^,,.}) - (T{vE',,,.^vl/+^^^,}) . {^+a,.<^,,^.,a..) ■ 

(018) 

By an analysis similar to the one in Section fill D[ we get 

>5'2,1;m(z; X, y) = 5';^^[^,(z; X, y) + 5^1:^,(2; x, y) + S^^^^.^iz; x, y) , (019) 

where the first term is asymptotically negligible with respect to the last two for large distances. 
The terms 5^j.^(z; x, y) and S'^]^.^j(z; x, y) correspond to contributions to the correlation func- 
tion coming from the infrared integration, whose computation requires, as in Sections IIII CI and 
IIIIDI the decomposition of the infrared field into the sum of quasi-particle fields indexed by 
w = ± and supported, in momentum space, around the two different Fermi points p*^. By the 
compact support properties of the infrared fields, in the terms contributing to S'^]^.^(z; x, y), 
the quasi-particle indeces corresponding to the fields located at z are the same, and will be 
denoted by oj^; similarly, the quasi-particle indeces corresponding to the fields located at x and 
y are the same, and will be denoted by uj^y Finally, 5^;^.^(z; x, y) is defined as the sum of 
all the contributions such that lOz — ^xy, while 52^]^.^(z; x, y) corresponds to the terms with 
Wz = —LJxy By construction, 5^]^.^(z; x, y) can be written as a sum over two terms: 

5+i,,,(z;x,y) = ^ e-^^^-^'-y^ S+,.^^Jz;K,y) 020) 

Ul=± 
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and we can combine such terms in a single matrix 




021) 

It is clear from the multiscale construction of these correlation functions that, with a proper 
choice of the parameters, such matrices are asymptotically close to the Schwinger function of 
the QED2+1 model seen above, as explained by the following theorem, which is, in fact, a 
corollary of the analysis in the previous sections and of a finite dimensional fixed point argument. 

Theorem 3. Let U and eo be small enough. It is possible to choose Z and c in i fCl J|) - fC15|j 
as functions ofU,eo,M and v{0), so that, i/k, p,k — p are all on the same scale h (i.e., if 
ao7 



h-l 



< |k|,|p|,|k-p| <ao7''+S /^<0;, 



e(k) = (Mk) (1 + 0(7''')) , 022) 

r,,(k,p) = V-k^k+p) (1 + 0(7"')) , 023) 

where Zq,Z\,Zi in I CI ^.?|] are suitable constants, depending on U,eo,M and v{0) and 
0<9<1. 



Theorem 3 says that, by choosing the wave function renormalization and the velocity of light 
in the QED model as suitable functions of U, eo, M and v{0), its two point Schwinger functions 
coincide with the ones of the Hubbard model, up to corrections which are negligible at small 
momenta. With this choice of Z and c, the vertex functions of QED2+1 are asymptotically 
proportional to those of the Hubbard model, provided that the renormalizations are properly 
chosen. Note that, while in a relativistic QFT is /i- independent, here it is not [ij], the 
symmetry (9') being broken by the underlying lattice; however, one can check, by arguments 
similar to the ones used in the proof of Lemma 2, that the lattice symmetries imply that the 
renormalizations Z^ are still diagonal in jj,: note, in fact, that in principle the r.h.s. of (|CI23|) 
could be of the form 

5]z^,.(j-;Vk^k+p)(l + 0(7"')), 

but, remarkably, Z^^^i, turns out to be diagonal. 

Theorem 3 implies that the Schwinger functions of the 2D Hubbard model on the honeycomb 
lattice obey to a Ward Identity analogous to (|CI12|1 . as it follows by combining (|CI12|1 with 
((Cr22)) - ((Cn3)) . see [11] ■ This is true not only in the free case f/ = (in which case the WI 
can be verified by a simple explicit computation) but also, remarkably, in the interacting case. 
Note that, with respect to the WI for QED, the WI for the Hubbard model is modified by 
the presence of some proportionality constants, which take into account both the relativistic 
renormalization of the charge and the fact that the Hubbard model breaks some relativistic 
symmetries. 

Let us conclude by remarking that, while here the WI can be proved a posteriori of the 
construction of the correlation functions, in the presence of Coulomb interactions the validity 
of an analogous WI is believed to play a crucial role in the construction of the model itself, as 
in one dimension p^ : in fact, in that case, the interparticle interaction becomes marginal in a 
RG sense [lH and the presence of Wis is a key ingredient in the control of the flow of the beta 
function equation, as in QED or in the Luttinger model. 
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